Magnetic coupling in CoCr 2 04 and MnCr 2 04: an LSDA+[7 study 
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We present a first principles LSDA+t/ study of the magnetic coupling constants in the spinel mag- 
nets CoCr2 04 and MnCr2C>4. Our calculated coupling constants highlight the possible importance 
of AA interactions in spinel systems with magnetic ions on both A and B sites. Furthermore, we 
show that a careful analysis of the dependence of the magnetic coupling constants on the LSDA+f/ 
parameters provides valuable insights in the underlying coupling mechanisms, and allows to obtain 
a quantitative estimate of the magnetic coupling constants. We discuss in detail the capabilities and 
possible pitfalls of the LSDA+[7 method in determining magnetic coupling constants in complex 
transition metal oxides. 



I. INTRODUCTION 

Geometrically frustrated spin systems exhibit a low- 
temperature behavior that is fundamentally different 
from conventional (non-frustrated) spin systems*^ The 
incompatibility between local interactions and global 
symmetry in geometrically frustrated magnets leads to 
a macroscopic degeneracy that prevents these systems 
from ordering. In some cases this degeneracy is lifted by 
further neighbor interactions or by a symmetry-breaking 
lattice distortion, resulting in ordered spin structures 
at temperatures that are significantly lower than what 
would be expected simply from the strength of the near- 
est neighbor interaction. Since usually several different 
ordered configurations with comparable energy exist in 
these systems, a very rich low temperature phase dia- 
gram can be observed. 

Recently, it has been found in various magnetic spinel 
systems (general chemical formula: AB2X4) that the geo- 
metrical frustration among the B sites in the spinel struc- 
ture can give rise to pronounced effects due to spin-lattice 
coupling. In ZnCr204 and CdCr204 the macroscopic de- 
generacy is lifted by a tetragonal lattice distortion, re- 
sulting in complicated non-collinear spin ordering]^ In 
addition, a pronounced splitting of certain phonon modes 
due to strong spin-phonon coupling has been found in 
ZnCr204^i^ Non-collinear spiral magnetic ordering at 
low temperatures has also been found in CoCr204 and 
MnCr204f£i£ where the presence of a second magnetic 
cation on the spinel A site lifts the macroscopic degener- 
acy. Such non-collinear spiral magnetic order can break 
spatial inversion symmetry and lead to the appearance 
of a small electric polarization and pronounced magneto- 
electric coupling&12- Indeed, dielectric anomalies at the 
magnetic transition temperatures have been found in 
polycrystalline CoCi^O^r^ an d recently a small electric 
polarization has been detected in single crystals of the 
same material^ Magnetic spinels therefore constitute a 
particularly interesting class of frustrated spin systems 
exhibiting various forms of coupling between their mag- 
netic and structural properties. Furthermore, both A 



and B sites in the spinel structure can be occupied by 
various magnetic ions and simultaneously the X anion 
can be varied between O, S, or Se. This compositional 
flexibility opens up the possibility to chemically tune the 
properties of these systems. 

To understand the underlying mechanisms of the var- 
ious forms of magneto-structural coupling, it is impor- 
tant to first understand the complex magnetic structures 
found in these systems. Such complex magnetic struc- 
tures can be studied using model Hamiltonians for inter- 
acting spin systems, which can be treated either clas- 
sically or fully quantum mechanically. For the cubic 
spinel systems, a theory of the ground state spin config- 
uration has been presented by Lyons, Kaplan, Dwight, 
and Menyuk (LKDM) 13 about 45 years ago. Using a 
model of classical Heisenberg spins and considering only 
BB and AB nearest neighbor interactions, LKDM could 
show that in this case the ground state magnetic struc- 
ture is determined by the parameter 
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which represents the relative strength between the two 
different nearest neighbor interactions Jbb an d Jab- 1A 
For u < u = 8/9 the collinear Neel configuration, i.e. all 
A-site spins parallel to each other and anti-parallel to the 
-B-site spins, is the stable ground state. For u > uq it was 
shown that a ferrimagnetic spiral configuration has the 
lowest energy out of a large set of possible spin configura- 
tions and that it is locally stable for uo < u < u" w 1.298. 
For u > u" this ferrimagnetic spiral configuration is un- 
stable. Therefore, it was suggested that the ferrimagnetic 
spiral is very likely the ground state for uq < u < u" , but 
can definitely not be the ground state for u > u" 

On the other hand it has been found that neutron scat- 
tering data for both CoCr204 and MnCr204 are well 
described by the ferrimagnetic spiral configuration sug- 
gested by LKDM, although a fit of the experimental data 
to the theoretical spin structure leads to values of u ss 2.0 
for CoCr204j£ and u w 1.6 for MnCr204, 7 which accord- 
ing to the LKDM theory correspond to the locally unsta- 
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ble regime. Surprisingly, the overall agreement of the fit 
is better for CoCr 2 04 than for MnCr 2 04, even though 
the value of u for CoCr204 is further within the unsta- 
ble region than in the case of MnCr2 04. From this it 
has been concluded that: i) the ferrimagnetic spiral is a 
good approximation of the true ground state structure 
even for u > u" , ii) that the importance of effects not 
included in the theory of LKDM is probably more sig- 
nificant in MnCr2C>4 than in CoCr204, and iii) that the 
ferrimagnetic spiral is indeed very likely to be the true 
ground state for systems with u < u < u"Z&12- 

Recently, Tomiyasu et al. fitted their neutron scat- 
tering data for CoC^CU and MnCr204 using a ferri- 
magnetic spiral structure similar to the one proposed by 
LKDM but with the cone angles of the individual mag- 
netic sublattices not restricted to the LKDM theory*^ 
As originally suggested by LKDM, they interpreted their 
results as indicative of a collinear Neel-like ferrimagnetic 
component exhibiting long-range order below Tq and a 
spiral component, which exhibits only short-range order 
even in the lowest temperature phase. 

In order to assess the validity of the LKDM theory 
and to facilitate a better comparison with experimental 
data, an independent determination of the magnetic cou- 
pling constants in these systems is very desirable. Den- 
sity functional theory (DFT, see Ref.ftq) provides an effi- 
cient way for the ab initio determination of such magnetic 
coupling constants that can then be used for an accurate 
modeling of the spin structure of a particular system. 
DFT also offers a straightforward way to investigate the 
effect of structural distortions on the magnetic coupling 
constants, and is therefore ideally suited to study the 
coupling between magnetism and structural properties. 

Traditionally, insulating magnetic oxides represent a 
great challenge for DFT-based methods due to the strong 
Coulomb interaction between the localized d electrons. 
However, recently the local spin density approximation 
plus Hubbard U (LSDA+C/) method has been very suc- 
cessful in correctly determining various properties of 
such strongly correlated magnetic insulators^ In par- 
ticular, it has been used for the calculation of mag- 
netic coupling constants in a variety of transition metal 
oxides^i^^^ 

Here we present an LSDA+C/ study of the magnetic 
coupling constants in the spinel systems CoCr 2 04 and 
MnCr204. The goal of the present paper is to provide ac- 
curate values for the relevant coupling constants in these 
two systems, in order to test the assumptions made by 
LKDM and to resolve the uncertainties in the interpreta- 
tion of the experimental data. In addition, we assess the 
general question of how accurate such magnetic coupling 
constants in complex oxides can be determined using the 
LSDA+C/ method. 

We find that in contrast to the assumptions of the 
LKDM theory, the coupling between the A site cations 
is not necessarily negligible, but that the general valid- 
ity of the LKDM theory should be better for CoCr2 04 
than for MnCr2 04, in agreement with what has been 



concluded from the experimental data. However, in con- 
trast to what follows from fitting the experimental data 
to the LKDM theory, the calculated u for CoCr204 is 
smaller than the value of u — 2.0 obtained from the ex- 
perimental fit. In addition, we show that by analyzing 
the dependence of the magnetic coupling constants on 
the LSDA+C/ parameters and on the lattice constant, 
the various mechanisms contributing to the magnetic in- 
teraction can be identified, and a quantitative estimate 
of the corresponding coupling constant can be obtained 
within certain limits. 

The present paper is organized as follows. In Sec. |TT] 
we present the methods we use for our calculations. In 
particular, we give a brief overview over the LSDA+C/ 
method and the challenges in using this method as a 
quantitative and predictive tool. In Sec. IIIII we present 
our results for the lattice parameters, electronic struc- 
ture, and magnetic coupling constants of the two investi- 
gated Cr spinels. Furthermore, we analyze in detail the 
dependence of the magnetic coupling constants on the 
lattice constant and LSDA+C/ parameters, and we dis- 
cuss the reasons for the observed trends. We end with a 
summary of our main conclusions. 

II. METHODS 

A. LSDA+C/ 

The LSDA+C/ method offers an efficient way to cal- 
culate the electronic and magnetic properties of complex 
transition metal oxides. The idea behind the LSDA+C/ 
method is to explicitly include the Coulomb interaction 
between strongly localized d or / electrons in the spirit 
of a mean-field Hubbard model, whereas the interactions 
between the less localized s and p electrons are treated 
within the standard local spin density approximation 
(LSDA). 22 To achieve this, a Hubbard-like interaction 
term Ejj , which depends on the occupation of the local- 
ized orbitals, is added to the LSDA total energy, and an 
additional double counting correction E dc is introduced 
to subtract that part of the electron-electron interaction 
between the localized orbitals that is already included in 
the LSDA: 

E = £Lsda + E[j — Edc ■ (2) 

Here 

Ejj = — y ' (£^71737274 — ^71737472) n 7l72 n 7374 (3) 
{7} 

and 

E dc = -n{n-l)-—J2ns(n s -l) , (4) 

where 7 = (to, s) is a combined orbital and spin in- 
dex of the correlated orbitals, the correspond- 
ing orbital occupancy matrix, n s = ^ m n rns ^ ms and 
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n = J2 S n s are the corresponding traces with respect to 
spin and both spin and orbital degrees of freedom, and 

737274 — (nii'm>3\Vee\m'2ni4)5 SlS2 Ss 3 s 4 are the matrix 
elements of the screened electron electron interaction, 
which are expressed as usual in terms of two parame- 
ters, the Hubbard U and the intra-atomic Hund's rule 
parameter J (see Ref. Il7l ) . 

The LSDA+C/ method has been shown to give the cor- 
rect ground states for many strongly correlated magnetic 
insulators, and thus represents a significant improvement 
over the LSDA for such systems^ Furthermore, the 
LSDA+C/ method is very attractive due to its simplicity 
and the negligible additional computational effort com- 
pared to a conventional LSDA calculation, ft therefore 
has become a widely used tool for the study of strongly 
correlated magnetic insulators. Since the LSDA+C/ 
method treats the interactions between the occupied or- 
bitals only in an effective mean-field way, it fails to de- 
scribe systems where dynamic fluctuations are impor- 
tant. For such systems, the local density approximation 
plus dynamical mean field theory (LDA+DMFT), which 
also includes local dynamic correlations, has been intro- 
duced recently^ However, the LDA+DMFT method is 
computationally rather demanding, and is currently too 
costly to be used for the calculation of magnetic charac- 
teristics of such complex materials as the spinels. On the 
other hand, for a large number of systems such fluctua- 
tions are only of minor importance, and for these systems 
the LSDA+C/ method leads to a good description of the 
electronic and magnetic properties. 

However, in order to obtain reliable results, the use 
of the LSDA+C/ method should be accompanied by a 
careful analysis of all the uncertainties inherent in this 
method. An additional goal of the present paper is 
therefore to critically assess the predictive capabilities of 
the LSDA+C/ method for the determination of magnetic 
coupling constants in complex transition metal oxides. 
Apart from the question about the general applicabil- 
ity of the LSDA+C/ approach to the investigated sys- 
tem, and the unavoidable ambiguities in the definition of 
the LSDA+C/ energy functional (Eqs. ©-g]))^^ the 
proper choice of the parameters U and J H represents 
one of the main hurdles when the LSDA+C/ method is 
used as a quantitative and predictive tool. 

U and J H can in principle be calculated using con- 
strained density functional theory,— thus rendering the 
LSDA+C/ method effectively parameter- free. In practice 
however, the exact definition of U and J H within a solid 
is not obvious, and the calculated values depend on the 
choice of orbitals or the details of the method used for 
their determination. 27 ' 28 ' 29 ' 30 Therefore, parameters ob- 
tained for a certain choice of orbitals are not necessarily 
accurate for calculations using a different set of orbitals. 

In the present work we thus pursue a different ap- 
proach. We choose values for U and J H based on a 
combination of previous constrained DFT calculations, 
experimental data, and physical reasoning, and these val- 
ues are then varied within reasonable limits to study the 



resulting effect on the physical properties. In particular, 
for the spinel systems studied in this work the Hubbard 
Us on the transition metal sites are varied between 2 eV 
and 6eV (in leV increments), with the additional re- 
quirement that C/cr < Ua {A — Co, Mn). For the on-site 
Hund's rule coupling we use two different values, J H 
= OeV and J H = I eV with jg = jf. The condi- 
tions C/cr < Ua and Jq t = are motivated by con- 
strained DFT calculations for a series of transition metal 
perovskite systems, which showed that the Hubbard U 
increases continuously from V to Cu, whereas the on-site 
exchange parameter J H is more or less constant across 
the series^ A similar trend for U can be observed in 
the simple transition metal monoxides i 22 ' 29 Although in 
the spinel structure the coordination and formal charge 
state of the A cation is different from the B cation, we as- 
sume that the assumption C/o < Ua is nevertheless valid, 
since the screening on the sixfold coordinated B site is 
expected to be more effective than on the tetrahedral A 
site. Further evidence for the validity of this assump- 
tion is given by the relative widths of the d bands on 
the A and B sites obtained from the calculated orbitally 
resolved densities of states (see Fig. Q] and Sec. MI Bj) . 

The absolute values of U used in this work are moti- 
vated by recent constrained DFT calculations using lin- 
ear response techniques ; 29 ' 30 which lead to significantly 
smaller values of U than previous calculations using the 
linear muffin tin orbital (LMTO) method, where the oc- 
cupation numbers are constrained by simply setting all 
transfer matrix elements out of the corresponding or- 
bitals to zero i 28 ' 31 Typical values obtained for various 
transition metal ions in different chemical environments 
are between 3-6 eVi 29 ' 30 For the Cr 3+ ion a value of 
U ~ 3eV, derived by comparing the calculated densi- 
ties of states with photo-emission data, has been used 
successfully.— ^ We thus consider the values Ua = 4-5 eV 
and C/cr = 3 eV as the most adequate U parameters for 
our systems. Nevertheless, we vary these parameters here 
over a much larger range, in order to see and discuss the 
resulting trends in the calculated magnetic coupling con- 
stants. 

For the Hund's rule parameter J H screening effects 
are less important, and calculated values for various sys- 
tems are all around or slightly lower than 1 eV . 22 ' 28 On 
the other hand, a simplified LSDA+C/ formalism is some- 
times used, where the only effect of J H is to reduce the 
effective Coulomb interaction U e g — U - j H 33^3^35 j n 
this work we use the two values J H = eV and J H = 
1 eV to study the resulting effect on the magnetic cou- 
pling constants. 



B. Other technical details 

To determine the magnetic coupling constants corre- 
sponding to the closest neighbor magnetic interactions 
between the various sublattices, we calculate the total 
energy differences for four different collinear magnetic 
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configurations: the Neel type ferrimagnetic order, the 
ferromagnetic configuration, and two different configura- 
tions with anti-parallel magnetic moments within the A 
and B sub-lattices respectively, and we then project the 
resulting total energies on a simple classical Heisenberg 
model, 

E = -^Ji j Si-Sj = -^J i je i -e.j , (5) 

where only the nearest neighbor coupling constants Jab , 
Jbb, and J aa are assumed to be nonzero, and where we 
defined the coupling constants Jij — JijSiSj correspond- 
ing to normalized spin directions e, of the magnetic ions. 
We note that even though for itinerant systems such as 
the elementary magnets Fe, Co, and Ni, the coupling 
constants obtained in this way can be different from the 
ones obtained for only small variations from the collincar 
configurations^ the local magnetic moments of many in- 
sulating transition metal oxides, in particular the systems 
investigated in the present study, behave much more like 
classical Heisenberg spins and thus the simpler approach 
pursued in this work is justified. We point out that a de- 
termination of all possible further neighbor interactions 
is beyond the scope of this paper and is therefore left for 
future studies. 

We perform calculations at both experimentally de- 
termined lattice constants and theoretical lattice param- 
eters. The theoretical lattice parameters are obtained 
by a full structural relaxation within the LSDA for a 
collinear Neel-type magnetic configuration. The same 
LSDA lattice parameters are used in all our calcula- 
tions with varying values of the LSDA+t/ parameters 
U and J H . In order to reduce the required computa- 
tional effort, we do not perform relaxations for each indi- 
vidual set of LSDA+J7 parameters. Except when noted 
otherwise, all calculations are performed using the "Vi- 
enna Ab-initio Simulation Package" (VASP) employing 
the projector augmented wave (PAW) method i 37 i 38 ' 39 We 
use a plane wave energy cutoff of 450 eV (550 eV for re- 
laxations) and a 5x5x5 T-centered mesh for Brillouin 
zone integrations. Increasing the mesh density by us- 
ing a 8x8x8 mesh results only in negligible changes for 
the calculated total energy differences. Structural re- 
laxations are performed until the forces are less than 
10~ 5 eV/A and all components of the stress tensor are 
smaller than 0.02 kbar. The electronic self-consistency 
cycle is iterated until the total energy is converged better 
than 10 -8 eV. In addition, we perform some test calcu- 
lations using the full-potential linear-augmented-plane- 
wave (FLAPW) method. 4° For these calculations we use 
the Wien97 code^i with our own implementation of the 
LSDA+f/ method. The plane-wave cut-off parameter is 
set to 223 eV in these calculation, and the Brillouin-zone 
integration is also carried out on a 5x5x5 T-centered 
mesh. The criterion for self-consistency is the difference 
in the total energy after the last two iterations being less 
than 10- 4 Ry. 



TABLE I: Structural parameters calculated in this work, a 
is the lattice constant of the cubic spinel structure, and the 
internal structural parameter x corresponds to the Wyckoff 
position 32e (x,x,x) of the oxygen sites. Columns "theo." con- 
tain the values calculated in this work while columns "exp." 
contain experimental data. 





CoCr 2 4 




MnCr 2 4 






exp. (Ref.Jl) 


theo. 


exp. (Ref.42) 


theo. 


a [A] 


8.335 


8.137 


8.435 


8.242 


X 


0.264 


0.260 


0.264 


0.262 



III. RESULTS AND DISCUSSION 

A. Structural relaxation 

Table [J shows the structural parameters obtained in 
this work together with corresponding experimental data. 
The theoretical lattice constants are obtained within the 
LSDA and for Neel-type ferrimagnetic order, and are 
about 2.3 % smaller than the corresponding experimen- 
tal values for both materials. The calculated internal 
structural parameters x are in very good agreement with 
experiment. The underestimation of the lattice constant 
by a few percent is a typical feature of the LSDA in com- 
plex transition metal oxides^^ 3 - 



B. Electronic structure 

Fig. [U shows the densities of states for both C0CT2O4 
and MnCr2 04 calculated using the LSDA and the 
LSDA+f/ method with U A = 5eV, U Cl = 3eV, J% = 
J§ T = eV, and a collinear Neel-type magnetic configura- 
tion at the experimental lattice constants. Both systems 
are insulating within the LSDA, although the LSDA en- 
ergy gap for CoCr2C>4 is very small, about 0.15 eV. The 
LSDA gap is larger for MnCr2 04, since in this system 
the gap is determined by the relatively strong crystal- 
field splitting on the octahedral B site and the equally 
strong magnetic splitting, whereas in C0CT2O4 the width 
of the LSDA gap is limited by the small crystal-field split- 
ting on the tetrahedrally coordinated A site. The use 
of the LSDA+t/ method increases the width of the en- 
ergy gap substantially and pushes the majority d states 
on the A site down in energy, leading to strong overlap 
with the oxygen 2p states. In the LSDA the transition 
metal d states are well separated from the oxygen p man- 
ifold, whereas the LSDA+C7 method increases the ener- 
getic overlap between these states. In all cases the gap 
is between occupied and unoccupied transition metal d 
states. 

It can be seen that the bandwidth of the <i-bands for 
the tetrahedrally coordinated A site is indeed smaller 
than for the octahedral B site. Thus, the d states on the 
A sites are more localized and one can expect a larger on- 
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FIG. 1: Densities of states (in states/eV) for CoCr204 (left 
two panels) and MnC^CU (right two panels) calculated us- 
ing the LSDA (upper two panels) and the LSDA+[/ method 
with Ua = 5 eV, U Cl = 3 eV, and jf = jg t = eV (lower 
two panels). Calculations were done at the experimental lat- 
tice parameters for a collinear Neel-type ferrimagnetic struc- 
ture where the direction of the Cr magnetic moments was de- 
fined as "spin-up" , and corresponding "spin-down" states are 
shown with a negative sign. The gray shaded areas represent 
the total density of states, the curves shaded with diagonal 
lines represent the d states on the A site of the spinel lattice, 
and the thick black lines correspond to the Cr d states. Zero 
energy separates the occupied from the unoccupied states. 
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FIG. 2: Magnetic coupling constants Jab (upper panels), 
Jbb (middle panels), and J a a (lower panels) calculated for 
CoCr2 04 (left) and MnCr2 04 (right) as a function of Ua for 
Ucr = 2eV (open circles), 3eV (filled squares), 4eV (open 
diamonds), and 5eV (filled triangles). All calculations were 
performed using experimental lattice parameters and J a = 
Jg = leV. 



site Coulomb interaction than on the Cr B site, in agree- 
ment with the assumption that Uci < Ua (see Sec. Ill At . 



C. Magnetic coupling constants 

Fig. [2] shows the magnetic coupling constants calcu- 
lated using the experimental lattice parameters, Sa = 
Jq t = 1 eV, and different values of the Hubbard U on 
the A and B sites. All coupling constants are nega- 
tive, i.e. antiferromagnetic, and decrease in strength 
when the Hubbard parameters are increased. The "inter- 
sublattice" coupling Jab depends similarly on both Ua 
and Ub, whereas both "intra-sublattice" coupling con- 
stants Jbb and Jaa depend only on the Hubbard pa- 
rameter of the corresponding sublattice. 

The BB interaction in the spinel lattice is known 
to result from a competition between antiferromagnetic 
(AFM) direct cation-cation exchange and indirect cation- 
anion-cation exchange, which for the present case of a 
90° cation-anion-cation bond angle gives rise to a fer- 
romagnetic (FM) interaction.— The AFM direct interac- 
tion is expected to dominate at smaller volumes, whereas 
at larger volumes the FM indirect interaction should be 
stronger. Furthermore, it is important to note that even 



the pure direct cation-cation interaction is comprised out 
of two parts: (i) the "potential exchange" due to the stan- 
dard Hcitlcr-London exchange-integral, which is always 
FM for orthogonal orbitals but is usually negligible, and 
(ii) the AFM "kinetic exchange" , which results from a 
second order perturbation treatment of the electron hop- 
ping and is proportional to 1/U , 44 i 45 The observed U 
dependence of Jbb can thus be understood as follows: 
At small values of U the AFM direct kinetic exchange 
is strongest, but it is suppressed as the value of U is in- 
creased. The FM indirect cation-anion-cation exchange 
also decreases, but in addition increasing U shifts the 
cation d states down in energy and thus leads to enhanced 
hybridization with the anion p states (see Sec. IIIIBI) . 
This enhanced hybridization partially compensates the 
effect of increasing U so that the indirect exchange de- 
creases slower than 1/U. Therefore, the FM indirect ex- 
change is less affected by increasing U than the AFM 
direct exchange, and thus gains in strength relative to 
the latter. This explains why the observed decrease of 
Jbb is stronger than 1/U. In fact, for the larger exper- 
imental volumes and using J = eV for the Hund's rule 
coupling (see discussion below) the BB coupling in both 
systems even becomes slightly FM for large U. 
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The AB coupling in the spinels is mediated by a cation- 
anion-cation bond with an intermediate angle of ~120°, 
which makes it difficult to predict the sign of the coupling 
based on general considerations. A weak AFM interac- 
tion has been proposed for the case of empty e g orbitals 
on the B site^ in agreement with the present results. 

Comparing the values of Jab and Jb b calculated for a 
constant set of LSDA+L 7 parameters shows that both 
are of the same order of magnitude. On the other 
hand Jaa is expected to be significantly weaker, since 
it corresponds to a cation-cation distance of ~ 3.6 A 
with the shortest superexchange path along an A-O-O- 
A bond sequence. Based on this assumption 3 aa was 
neglected in the theoretical treatment of LKDM. 13 Nev- 
ertheless, in our calculations J aa is found to be of appre- 
ciable strength. This is particular striking for MnCr204, 
but also for CoC^CU the difference between Jaa and 
Jab {Jbb) is less than an order of magnitude. From 
this it follows that Jaa is definitely non-negligible in 
MnCr204 and can also lead to significant deviations from 
the LKDM theory in CoCr204. We point out that this 
conclusion holds true independently of the precise val- 
ues of the LSDA+[/ parameters used in the calculation. 
An appreciable value for Jaa has also been found in a 
previous LSD A study of the spinel ferrite MnFe2044£ 

As stated in Sec. IIIB1 a full determination of all pos- 
sible further neighbor interactions is beyond the scope 
of this paper. However, to obtain a rough estimate of 
the strength of further neighbor interactions in CoCr204, 
and to see whether this affects the values of Jab, Jbb, 
and Jab obtained in this work, we perform some addi- 
tional calculations using a doubled unit cell. This al- 

(3) 

lows us to determine the coupling constant J BB , corre- 
sponding to the third nearest neighbor BB coupling. As 
shown in Ref.H^, due to the special geometry of the spinel 
structure, this third nearest neighbor coupling is larger 
than all other further neighbor interactions within the 
B sublattice, and can be expected to represent the next 
strongest magnetic interaction apart from Jbb, Jab, and 
Jaa ■ This coupling is mediated by a B-O-O-B bond se- 
quence and corresponds to a BB distance of 5.89 A. For 
comparison, the distances corresponding to Jbb, Jab, 
and Jaa are 2.94 A, 3.46 A, and 3.61 A, respectively^ 
For these test calculations we use the experimental lat- 
tice parameters of CoCr204 and the LSDA+L 7 parame- 
ters U Co = 5 eV, U G T = 3 eV, and J H = eV. We obtain 

(3) 

a value of J BB — 0.15 eV, corresponding to a weak FM 

(3) 

coupling. However, the magnitude of J BB is small com- 
pared to Jab and Jbb, and we therefore continue to 
neglect further neigbor interactions in the following. 

Next we calculate the magnetic coupling constants 
using the lattice parameters obtained by a full struc- 
tural optimization within the LSDA, and also by using 
J H = eV at both experimental and theoretical lattice 
parameters. Again, we vary Ua and Uct independently. 
The observed variation of the coupling constants with re- 
spect to the Hubbard parameters is very similar to the 



TABLE II: Calculated magnetic coupling constants Jab, 
Jbb, and Jaa for different lattice parameters and different 
values of the intra-atomic Hund's rule coupling parameter J H 
for Ua = 5 eV and Ucr = 3 eV. Lattice parameters "exp." and 
"theo." refer to the corresponding values listed in Table U 





J H (eV) 


1.0 


0.0 


1.0 


0.0 




latt. param. 


exp. 


exp. 


theo. 


theo. 




Jab (meV) 


-4.44 


-3.55 


-6.02 


-4.83 


CoCr 2 4 


Jbb (meV) 


-3.33 


-1.04 


-6.90 


-4.34 




Jaa (meV) 


-0.50 


-0.44 


-0.77 


-0.58 




Jab (meV) 


-3.14 


-1.40 


-4.88 


-2.61 


MnCr 2 4 


Jbb (meV) 


-2.91 


-0.74 


-5.22 


-2.74 




Jaa (meV) 


-1.19 


-0.92 


-1.88 


-1.45 



case shown in Fig. [2j only the overall magnitude of the 
magnetic coupling constants is changed. We therefore 
discuss only the results obtained for Ucr = 3eV and 
Ua — 5eV, which are physically reasonable choices for 
these parameters, as discussed in Sec. Ill Al 

Table [II] shows the calculated magnetic coupling con- 
stants for the different cases. It is apparent that both vol- 
ume and the intra-atomic exchange parameter J H have 
a significant effect on the calculated results. The vol- 
ume dependence can easily be understood. The smaller 
theoretical volume leads to stronger coupling between 
the magnetic ions. This is particularly significant for 
Jbb, since the direct exchange interaction between the 
B cations is especially sensitive to the inter-cation dis- 
tance. The corresponding coupling is therefore strongly 
enhanced (suppressed) by decreasing (increasing) the lat- 
tice constant. The indirect superexchange interaction 
also depends strongly on the inter-atomic distances. 

It can be seen from Table [H] that J H — significantly 
decreases the strength of all magnetic coupling constants 
compared to J H = 1 eV. A strong J H dependence of 
the magnetic coupling has also been observed in other 
Cr spinels with non-magnetic cations on the A site^ 
Further calculations, with different values for J H on the 
A and B sites respectively, show that it is mostly Jq t 
which is responsible for this effect. On the other hand, 
varying Jj± nas a smaller effect on the magnetic coupling. 
This is consistent with the very strong J H dependence 
of Jbb and the weaker J H dependence of Jaa seen in 
Table El 

To understand the strong effect of jff. on the magnetic 
coupling constants we first take a look at the occupation 
numbers of the Cr d orbitals. The correspond- 

ing occupation numbers in CoCr204 are (calculated for 
a FM configuration at the experimental lattice parame- 
ters and using J H = 0): «t 2g ,| = 0.95, n t2g ,i = 0.05, 
n e g ,] — 0.32, and n eg ,i = 0.21. As expected, the oc- 
cupation of the t2 g orbitals represents the formal d va- 
lency with full spin-polarization, but in addition there is 
a sizable e g occupation, which contributes ~ 0.2/xb to 
the local spin moment of the Cr cation. This partial e g 
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occupation, which is due to hybridization with the oxy- 
gen p bands, gives rise to a FM interaction between the 
Cr sites, because the e g polarization is coupled to the 
t2 g spins via the Hund's rule coupling. 44 This FM inter- 
action between the Cr sites should therefore be propor- 
tional to the strength of the Hund's rule coupling. Thus, 
the stronger AFM interaction for J H = leV compared 
to J H = (see Table H3J might be surprising at first. 
However, it is important to realize that even though the 
parameter J H represents the strength of the Hund's rule 
coupling, its effect within the LSDA+17 framework is not 
to introduce a strong Hund's rule interaction. If one ana- 
lyzes the LSDA+f/ energy expression, Eq. ([2]), in a some- 
what simplified picture where the occupation matrix is 
diagonal and the Coulomb matrix elements arc orbitally 
independent, one can see that the double counting cor- 
rection, E'dc, exactly cancels the different potential shifts 
for orbitals with parallel and antiparallel spins that are 
caused by E v for J H ^ 0, if one of the d orbitals is filled. 
Thus, Eu — Ea c does not lead to an additional Hund's 
rule interaction compared to -Elsda, even for J ^ 0. 
It is generally assumed that this type of interaction is al- 
ready well described on the LSDA level. The only effect 
of J H is therefore an effective reduction of the on-site 
Coulomb repulsion. This can be seen in the following, 
where we write the simplified version of Eq. ([2|) as (see 
Ref.il): 

E = -Elsda + -^y- - n^n^j , (6) 

with U cfi = U - J H . Within this simplified LSDA+J7 
version, the effect of J H on the magnetic coupling con- 
stant can be understood as an effective reduction of the 
on-site Coulomb interaction. According to the previ- 
ously discussed U dependence of the magnetic coupling 
constants (see also Fig. [2]), a reduced on-site Coulomb 
interaction leads to a stronger AFM interaction for all 
calculated magnetic coupling constants. 

From Fig. [2] it can be seen that the magnetic cou- 
pling constants for CoCr02 using experimental lattice 
parameters and the LSDA+C/ parameters Uc = 6eV, 
U C t = 4eV, and J H = I eV, i.e. U cS ,co = 5eV and 
f7 effCr = 3eV, are: J AB = -3.26eV, J BB = -2.12 eV, 
and Jaa = — 0.30eV. The corresponding result for Uq — 
5eV, Uqi — 3eV, and J H = 0, i.e. for the same values 
of U e ff but different J H , are: Jab = — 3.55eV, Jbb = 
-1.04 eV, and J A a = -0.44eV (see Table QTJ. Thus, the 
pure dependence on U e ff seems to be approximately valid 
for Jab and Jaa, whereas there is a notable quantita- 
tive deviation from the simplified LSDA+[/ model in the 
case of Jbb- Nevertheless, the overall trend can still be 
understood from the simplified LSDA+t/ pictue. 

Finally, to assess the possible influence of different 
methods to solve the self-consistent Kohn-Sham equa- 
tions on the calculated magnetic coupling constants, we 
perform additional tests using a different electronic struc- 
ture code employing the FLAPW method (see Sec. Ill B| ). 



TABLE III: Magnetic coupling constants of CoC^CU and 
MnCr2 04 calculated using two different methods (FLAPW 
and PAW), different values for J , and the experimental lat- 
tice parameters. 







FLAPW 


PAW 






J H (eV) 


0.0 


1.0 


0.0 


1.0 




Jab (meV) 


-3.62 


-4.32 


-3.55 - 


•4.44 


CoCr 2 4 


Jbb (meV) 


-1.32 


-3.09 


-1.04 - 


3.33 




Jaa (meV) 


-0.23 


-0.00 


-0.44 - 


-0.50 




Jab (meV) 


-1.73 


-3.23 


-1.40 - 


-3.14 


MnCr 2 4 


Jbb (meV) 


-1.32 


-3.21 


-0.74 - 


-2.91 




Jaa (meV) 


-0.67 


-1.06 


-0.92 - 


1.19 



The results are summarized and compared to the PAW 
results in Table Hill There are some variations in the ab- 
solute values of the magnetic coupling constants obtained 
with the two different methods, but overall the agreement 
is rather good. Trends are the same in both methods, and 
in particular the strong effect of the LSDA+J7 Hund's 
rule parameter J H on the magnetic coupling constants 
is confirmed by the FLAPW calculations. One possible 
reason for the differences between the PAW results and 
the FLAPW results is that the radii of the projection 
spheres used in the PAW method are chosen differently 
from the radii of the Muffin-Tin spheres used to construct 
the FLAPW basis functions. 



D. The LKDM parameter u 

Figure [3] shows the variation of the LKDM parameter 
u = 4JbbSb = ^ Jbb with the strength of the on-site 

3JabSa 3 Jab _ b 

Coulomb interactions for the different lattice parameters 
and values of J H used in this work. The behavior of 
u follows from the corresponding trends in the coupling 
constants Jab and Jb b discussed in the previous section. 
Increasing Ua decreases the strength of Jab but leaves 
Jbb more or less unchanged, and thus increases the value 
of u. On the other hand both Jab and Jbb decrease 
with increasing Uqt, but the decrease is stronger for Jbb 
and therefore u decreases with increasing Uc T . Thus, the 
trends caused by the Hubbard parameters corresponding 
to the two different magnetic sites are opposite to each 
other. 

As already pointed out in the previous section, chang- 
ing the value of the LSDA+Z7 parameter J H and using 
different lattice constants essentially just shifts the over- 
all scale for the magnetic coupling constants without al- 
tering their U dependence. Therefore, using the larger 
experimental volume decreases u compared to the value 
obtained at the theoretical lattice parameters due to the 
very strong volume dependence of Jbb- Introducing the 
on-site Hund's rule coupling J H increases u, since Jbb 
is stronger affected by this and thus increases relative to 
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FIG. 3: Dependence of the LKDM parameter u on the Hub- 
bard U parameters of the different magnetic cations. Left 
panels correspond to CoC^CU, right panels to MnC^CU. 
From top to bottom the different panels correspond to cal- 
culations for exp. volume and J = 1 eV, exp. volume and 
J H = OeV, theo. volume and J H — leV, and theo. vol- 
ume and J H = OeV (open circles: Ucr = 2eV, filled squares: 
Ucr = 3 eV, open diamonds: Ucr = 4 eV, filled triangles: 
Uci = 5 eV) . Dashed horizontal lines indicate the critical val- 
ues u = 8/9 and u" w 1.298. 



Jab- 

For similar values of U a and Uqi the LKDM parame- 
ter u is larger in MnCr204 than in CoCr204, except for 
J H = 1 eV at the theoretical lattice parameters, where 
they are approximately equal. This is in contrast to what 
has been concluded by fitting the experimental neutron 
spectra to the spiral spin structure of the LKDM the- 
ory, which leads to the values u=1.6 for MnC^C^ and 
u=2.0 for CoCr 2 4 ,-£ i.e. the fitted value for CoCr 2 4 
is significantly larger than the value for MnCr204. 

We now try to give a quantitative estimate of u in 
the two systems. The first question is whether using ex- 
perimental or theoretical lattice constants leads to more 
realistic magnetic coupling constants. This question is 
not easy to answer in general. On the one hand, the 
LSDA underestimation of the lattice constant can lead 
to an overestimation of the magnetic coupling, since the 
cations are too close together and can therefore inter- 
act stronger than at the experimental volume. On the 
other hand, the indirect cation-anion-cation interaction 



is intimately connected to the chemical bonding^ 4 If the 
larger experimental lattice constant is used, this bonding 
is artificially suppressed and the corresponding magnetic 
coupling is eventually underestimated. It is therefore not 
obvious whether it is better to calculate the coupling con- 
stants at the experimental or the theoretical lattice pa- 
rameters, but the two cases at least provide reasonable 
limits for the magnetic coupling constants. We note that 
using the LSDA+t/ method for the structural relaxation, 
usually leads to lattice parameters that are in slightly 
better agreement with the experimental values, 43 which 
will decrease the corresponding uncertainty in the mag- 
netic coupling constants. In the present paper we do not 
perform such structural relaxations for each combination 
of the LSDA+t/ parameters, in order to reduce the re- 
quired computational effort. In addition, this allows us 
to discuss the pure effect of U and J H on the magnetic 
coupling constants, without contributions due to varying 
lattice parameters. 

Fig. [3] shows that for the physically reasonable param- 
eters U C r = 3 eV, U A = 4-5 eV, and J H = 1 eV the value 
of u in CoCr204 calculated at the theoretical lattice con- 
stant is slightly larger than the critical value u" « 1.298, 
where within the LKDM theory the ferrimagnetic spiral 
configuration becomes unstable. In MnCi'204 the cor- 
responding value is about equal to u" . At the experi- 
mental lattice constants the values of u in both systems 
are smaller than at the theoretical lattice constants, with 
the stronger effect in CoCr204, where u at the theoreti- 
cal lattice constant is about equal to uq = 8/9, the value 
below which, according to LKDM, a collinear ferrimag- 
netic spin configuration is the ground state. In MnCi'204 
the value of u at the experimental lattice constant is be- 
tween Mo and u" . Thus, in all cases the calculated values 
of u are consistent with the experimental evidence for 
noncollinear ordering. 

Since in MnCi'20 4 the calculation predicts a rather 
strong Jaa, the validity of the LKDM theory is question- 
able for this system, but for CoCr20 4 , where the magni- 
tude of Jaa is indeed significantly smaller than both Jab 
and Jbb, this theory should at least be approximately 
correct. However, for CoCr204 the calculated u both 
at experimental and at the theoretical lattice constant 
(and using physically reasonable values for the LSDA+C7 
parameters) is still significantly smaller than the value 
u = 2.0 obtained by fitting the experimental data to 
the LKDM theory^ It would therefore be interesting to 
study how the incorporation of Jaa into a generalized 
LKDM theory alters the conclusions drawn from the ex- 
perimental data. Obviously, a non-negligible Jaa will 
further destabilize the collinear Neel configuration, but 
the possible influence of Jaa on the ferrimagnetic spi- 
ral structure cannot be obtained straightforwardly. Of 
course it cannot be fully excluded that the discrepancy 
between the calculated value of u for CoCr204 and the 
value extracted from the experimental data is caused by 
some deficiencies of the LSDA+C/ method. For exam- 
ple, it was shown in Ref. [H that for MnO the LSDA+J7 



9 



method does not offer enough degrees of freedom to cor- 
rectly reproduce both nearest and next nearest neighbor 
magnetic coupling constants. 

Finally, we note that the fact that J aa is not negli- 
gible in MnCr2 04 but has a significantly smaller magni- 
tude than Jab and Jbb in CoCr204 is compatible with 
the fact that the overall agreement between the experi- 
mental data and the LKDM theory is better for CoCr2 04 
than for MnC^C^. 8 However, a quantitative discrepancy 
between the value of u for CoCr 2 04 calculated in this 
work and the value derived from the experimental data 
remains. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented a detailed LSDA+f/ 
study of magnetic coupling constants in the spinel sys- 
tems CoCr204 and MnC^CU. We have found that the 
coupling between the A site cations, which is neglected 
in the classical theory of LKDM, is of appreciable size in 
CoCr2C>4 and definitely not negligible in MnCi^C^. The 
calculated LKDM parameter u, which describes the rel- 
ative strength of the BB coupling compared to the AB 
coupling and determines the nature of the ground state 
spin configuration in the LKDM theory, is found to be 
smaller than the values obtained by fitting experimental 
neutron data to the predictions of the LKDM theory. It 
remains to be seen whether this discrepancy is caused by 
the simplifications made in the LKDM theory, or whether 



it is due to deficiencies of the LSDA+C/ method used in 
our calculations. 

In addition, we have shown that it is difficult, but pos- 
sible, to arrive at quantitative predictions of magnetic 
coupling constants using the LSDA+Z7 method. In addi- 
tion, by analyzing the U and J H dependence of the mag- 
netic coupling constants it is possible to identify the var- 
ious interaction mechanisms contributing to the overall 
magnetic coupling. The presence of two different mag- 
netic cations with different charge states and different 
anion coordination, promotes the systems investigated 
in this work to a very hard test case for the predic- 
tive capabilities of the LSDA+t/ method. Nevertheless, 
some insight can be gained by a careful analysis of all 
methodological uncertainties, and the magnitudes of the 
magnetic coupling constants can be determined to a de- 
gree of accuracy that allows to establish important trends 
and predict the correct order of magnitude for the corre- 
sponding effects. 
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